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Abstract. In nonrelativistic quantum mechanics, the total (i.e. orbital plus 
spin) angular momentum of a charged particle with spin that moves in a Coulomb 
plus spin-orbit-coupling potential is conserved. In a classical nonrelativistic 
treatment of this problem, in which the Lagrange equations determine the orbital 
motion and the Thomas equation yields the rate of change of the spin, the 
particle's total angular momentum in which the orbital angular momentum is 
defined in terms of the kinetic momentum is generally not conserved. However, 
a generalized total angular momentum, in which the orbital part is defined in 
terms of the canonical momentum, is conserved. This illustrates the fact that 
the quantum-mechanical operator of momentum corresponds to the canonical 
momentum of classical mechanics. 



1. Introduction 

Quantum-mechanical analysis of atomic spectra is based on the conservation of the 
total angular momentum, defined as the sum of the orbital angular momentum and 
intrinsic spin, of an electron moving in a central plus the associated spin-orbit-coupling 
potential (see, e.g., [I]). This conservation follows immediately from the fact that, in 
nonrelativistic quantum mechanics, the operators of the Cartesian components and 
the square of the total angular momentum commute with the operator of the spin- 
orbit coupling and, thus, with the Hamiltonian of the system. However, the orbital 
angular momentum and the spin themselves are not conserved separately since the 
operators of their components do not commute with the spin-orbit-coupling operator 
(though the operators of their squares do). 

In this paper, we consider a classical nonrelativistic treatment of this problem 
using a Lagrangian in which the spin and orbital degrees of freedom of a charged 
particle with spin moving in a central electric field are coupled through a spin-orbit 
interaction Vso{r)s-l, where s and I — rx rav are the particle's spin and orbital angular 
momenta. The resulting Lagrange equations determine both the orbital motion and 
the motion of the spin. The equation of spin motion obtained is the nonrelativistic 
limit of the Thomas equation, derived by Thomas in 1927 in his historic, still essentially 
classical, analysis of the motion of a spinning electron in a Coulomb field [2] . In light 
of the quantum-mechanical result, it would not be unreasonable to expect that this 
classical treatment will yield the conservation of the total angular momentum I -\- s. 

It turns out, however, that, in this problem, the angular momentum Z -I- s is in 
general not conserved. But we find that the total angular momentum is conserved 
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in general when its orbital part is defined as L — r x P, where P is the canonical 
(generalized) momentum. When it is realized that canonical momentum is the classical 
counterpart of the quantum-mechanical operator of momentum, the classical and 
quantum-mechanical treatments are in harmony after all. 

The role of canonical momentum in the conservation laws of classical mechanics 
has been highlighted in a pedagogical paper earlier [3] , but a comparison with quantum 
mechanics does not seem to have been used in this context yet. Students typically 
encounter spin-orbit interaction only in courses on atomic physics and quantum 
mechanics, and thus an impression is easily created that such an interaction is a 
purely quantum phenomenon. A classical exposition of the spin-orbit interaction and 
its effect on the conservation of angular momentum should therefore be instructive. 
The present topic may serve as an illustration of the relationship between the results of 
classical and quantum-mechanical treatments of problems in an undergraduate course 
on atomic physics, assuming only some knowledge of the basic Lagrangian formalism. 



2. Classical equations of motion of a charged particle with spin in a 
Coulomb field 

Classical relativistic theory of combined spin and orbital motion (see, e.g., [4j [5l [6l 
El [HI 12]) is a difhcult subject in which one has to deal with complications such as 
that the centres of mass and charge of a relativistic spinning point particle do not 
coincide, which can lead to the so-called zitterbewegung ('jittery', or 'quivery' motion), 
identified in force- free solutions of the Dirac equation by Schrodinger [TD] , but observed 
experimentally only very recently, in a quantum simulation of the Dirac equation . 

However, a simple non-relativistic limit of a fully relativistic Lagrangian for 
a classical charged particle with spin can be derived. In this limit, the spin is 
parameterized in terms of vectors a and b as s = a x b, and the following Lagrangian 
is obtained for a classical particle of mass m, charge q and nonzero spin, moving with 
a velocity v in the electromagnetic field E, B with potentials A [7]: 



C=- mv'^ — q(f)+- v-A + - {b ■ d - b ■ a) + (axb) 
2 c 2 



qq / „ V \ vxv 
' B - -xE] - 



2mc V c / 2c^ 



■(1) 



Here, the dots on the symbols denote time derivatives and g is the so-called 
gyromagnetic factor ((;- factor for short), in terms of which the magnetic dipole moment 
/i of the particle is expressed as 

While both the magnitude of the spin s and the ^-factor are arbitrary, these quantities 
remain constant for a given classical particle; unlike in quantum mechanics, all the 
Cartesian components of the vectors s and /x are well defined. In order that Lagrangian 
([T]) can be used for setting up equations of motion, the acceleration v in its spin- 
dependent part is understood to be that which would arise in the absence of a spin- 
dependent interaction, i.e. v = {q/m){E + v x B/c) (see the discussion leading to 
equations (11.163) and (11.168) in [12]). 

Applying now Lagrangian ([1]) to a particle of charge g = — e (e > 0), moving 
in a Coulomb potential cj) = e/r, so that A = 0, E = — V(e/r) = er/r^ and 
B = 0, and taking accordingly v in the spin-dependent part of the interaction as 
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i) = —eE/m = —e r/mr gives 

1 1 ■ (q — l)e^ 

- mv^ H h-(b-d-b-a)- — — ^rr- (a x b) ■ (r x v) 

2 r 2 2mc^r'' 

1 , 1 



^.-^ . ^ . ^{b-d-b-a)-Vso{r)s-l. (3) 

Here, I = r x mv is the particle's orbital angular momentum defined in terms of the 
kinetic momentum mv, and Vso{r), the strength of the spin-orbit potential, is given 

by 

y^oir) = (4) 

For (7 = 2, which is to a very good approximation the factor of the electron, the 
spin-orbit strength has the empirically correct value that accounts for the fine 
structure of the spectrum of the hydrogen atom (see, e.g., [IB], section 34). 

Before we shall proceed with deriving the equations of motion from Lagrangian 
([3]|, a digression concerning the spin-orbit-interaction strength (|4]) is in order. 
According to classical electrodynamics, the coupling of the spin and orbital degrees of 
freedom arises through the interaction 

Ud = -d-E (5) 

of an electric field E with the electric dipole moment 

d=-(vxfi) (6) 

c 

that is relativistically acquired by a moving magnetic dipole /x [14l [Tsj^ . Putting 
E = er/r^ and using ([2]) to express the magnetic moment /x in terms of the spin s, 
interaction ([5]) is evaluated as 

Ud = — (f X /x) • — 

c 

= — ^/x- (r X -u) 

geP" 

However, for g — 2, the strength of this spin-orbit interaction is by a factor of 2 greater 
than the empirically correct value (|31). 

This discrepancy was removed in the early days of quantum mechanics by 
the famous Thomas precession factor of ^ [H [TB], which is an effect of the 
relativistic kinematics of curvilinear motion. Due to the noncommutativity of Lorentz 
transformations, two successive Lorentz boosts are in general equivalent to a single 
one plus a three-dimensional rotation. As a result, the rest frame of a particle moving 
along a curved path rotates with respect to the laboratory frame with an angular 
frequency that is given to first order in w/c as 

V X V 

-T-^. (8) 

t Interestingly, the Schwinger scattering 1161 . which is the scattering of neutrons by the electric 
field of an atomic nucleus, also arises through electric dipole moment ^ via interaction l[5j. The 
Schwinger-scattering Hamiltonian (see |13l . equation (42.1)) is transcribed classically in terms of the 
canonical momentum P as H = /2m — P ■ {fj, X E)/mc, which can be shown easily to be derivable 
from a Lagrangian L = mv^ /2 + {v X fi) • E/c | 17| . 
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(An instructive exposition of the Thomas precession can be found in |12) . section 
11.8.) In Lagrangian ([Ij, the Thomas precession is accounted for by the interaction 
term —{axb)-{vxv)/ 2(? = — s • Wt ■ This gives rise to the term — 1 in the factor g — \ 
of the empirically correct spin-orbit strength Q. As we shall see, all spin-dependent 
quantities will carry a factor (g — 1) instead of just g when the effect of the Thomas 
precession is included. 

We now turn to the particle's dynamics. With Lagrangian ([3|), we have 

dC, 

— = mvi + mVsoir){r x s)i, (9) 
ovi 

-— = — -Xi - m — r -{v X s) - mVso(r){v x sj^. (10) 

OXi r-^ dr r 

Here, the evaluation of the partial derivatives was facilitated by the identities 

s ■ {r X v) = V ■ {s X r) = r ■ {v X s). (11) 
The Lagrange equations 

d_dC_ _d£ j _ 1 2 3 fl2) 

At dvi dxi ' 

then give the following force on the particle, understood as mass times acceleration: 

e2 ^ dVso(r) „ d 

mv — — - f — m — ^ — [r • (v X s)]f — mVso{r){v x s) + —[rnVsoir^is x r)]. (13) 
dr at 

It is noteworthy that except for the Thomas-precession term —1 in the factor 

{g — 1) of the spin-orbit strength Vso{r)^ the spin-dependent part of force (fT3|) can 

be obtained by transforming the force on the particle's magnetic dipole from its 

instantaneous rest frame to the laboratory frame. Denoting instantaneous-rest-frame 

quantities by primes, the force on a magnetic dipole is given by |19l 120] 



= (M'-V')B'-ii^xf;' + — m'x J'. (14) 

c at' c 

where B' is the magnetic field and (/x' x E')lc is the so-called hidden momentum of 
the dipole (see, e.g., [3T] and references therein); the second line gives an equivalent 
expression obtained using a vector identity and the Ampere-Maxwell la'vili- In our 
case, the external electric current density J' vanishes at the dipole's location and, to 
first order in ti/c, i?' = — (t> x E)/c and E' = E, since the laboratory-frame magnetic 
field B is zercfjl; also, to first order in v/c, ti(/x'- V') = w(/x- V) and d/x'/d<' — d/x/d^^. 
Equation ^1)1 thus yields 

F' ^F + 0{v^/c^), (15) 



where 



-i^.V)[l.xE)+lEx'^. (16) 



§ Dropping the primes, V(/i-B) = (/x- V)B+/i X (V X B) (/x is not a function of spatial coordinates); 
V X B = AnJ/c + dE/cdt. 
II See [12], equation (11.149). 

^ This can be shown using equations (11.22) and (11.76) of |12) and equation (7) of |15) . (Note the 
different use of primes in the latter reference.) 
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(19) 



For the Coulomb field E = er/r^, the force can now be shown by straightforward 
vector algebra (see appendix A) to equal the spin-dependent part of the force ([T^ with 
Vso{r) = ge^ / (2m?c^r^) and the spin s expressed in terms of the magnetic moment /i 
using equation We note that the force on a moving magnetic dipole derived in 
[22] is given by the same expression as that of equation (IT6l) . 

The partial derivatives of Lagrangian ([3]) with respect to the spin variables and 
their time derivatives are 

1 ,• / N , . dC I , 

- = --b-Ko(r)6x^, - = -b (17) 

and 

The Lagrange equations for the spin degrees of freedom 

d_dC_dC d_dC_dC 

di'dd.^'da' di'db 'db 
therefore yield 

6 = Vsoir) Ixb, d = Vso{r) I x a, (20) 
and thus, since s = axb + axb, 

s = Vso{r) [{I X a) X b + a X (l X b)] 
= Vsoir) I X (axb) 

— mVso{r) (r X v) X s. (21) 

The rate of change of spin ([?!]) is that of the nonrelativistic limit of the Thomas 
equation [2] , which is a form of the well-known relativistic Bargmann-Michel-Telgedi 
(BMT) equation for an arbitrary classical spin and f;-factor [53] (instructive derivations 
of the BMT and Thomas equations are given in [T^], section 11.11). The Thomas 
equation to first order in v/c reads 

^'^ I B - -V X E] + — s X . (22) 



2mc \ c J mc 2c 

In our case, q = — e, E = er/r^ and B — 0. With these values, equation (15^ simplifies 
to 

s^ ^l'^Jf {rxv)xs, (23) 

which is indeed the same as equation (|21|) . recalling that the spin-orbit interaction 
strength Vso{r) in (f2T]) is given by (j!]). 

The physical content of the Thomas equation can be made more transparent 
by replacing qE in the second term of (|22p with mi), which is consistent with the 
nonrelativistic approximation employed, and then re-writing the equation in terms of 
magnetic moment ([2]) and the Thomas precession frequency ([8]): 

s = (1 X --VX E^ +u}T X s. (24) 

We recognize in the first term on the right-hand side the torque ^ x B' on a magnetic 
dipole fj, by its instantaneous-rest-frame magnetic field B' = B — {v x E)/c; the 
second term describes the Thomas precession of the spin. 
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3. Angular momentum of a classical charged particle with spin in a 
Coulomb field 

The rate of change of the orbital angular momentum defined in terms of the kinetic 
momentum mv as Z = r x mv equals the torque rxmv due to force (|13p : 

I ~ rxmv 

— —mVso{r)r x (v x s) + r x —\mVso(r)(s x r)] 

dt 

dVso(r) 

— —2mVso{r)r x {v x s) ~ ni — ^ — r x (r x s) — mVso{r)r x (r x s). (25) 

Here, the fact that r x f = was used. Employing explicit expression Q for Vso{r), 
we have 

3(g-l)e^ 

^ O 2 3 5 (26) 

the use of which in (|25|) yields 
(a - l)e^ 

I = -— ^ [2r X {v X s) - (3/r2)(r ■ v)r x {r x s) + r x {r x s)]. (27) 

Using now rates of change (1^51) and P7)) . it can be seen easily that 1 + s vanishes, 
i.e. the total angular momentum Z + s is conserved, only in the special case of the 
direction of the spin s being normal to the orbital plane defined by the vectors r and 
V, and, at the same time, the orbit being circular, i.e. r • v = 0. These are in fact the 
conditions for separate conservation of I and s. 

According to Noether's theorem, which associates a conserved quantity with each 
symmetry of the Lagrangian of a system (see, e.g., [21] )j there is a conserved angular 
momentum quantity if the system's Lagrangian is invariant with respect to a rotation. 
Lagrangian ([3]) is invariant under a rotation about the center of the Coulomb field; 
this fact can be shown to determine the conserved angular momentum quantity as 
L + s, where L is an orbital angular momentum defined in terms of the canonical 
momentum P = dC/dv as L = r x P. However, we here prefer to demonstrate the 
conservation of the angular momentum i + s by an explicit direct calculation, leaving 
the application of Noether's theorem to appendix B. 

Lagrangian ^ yields a canonical momentum 

dv 

^mv+ ^^ (rxs), (28) 

and the canonical orbital momentum is thus 
L = r X P 

(g - l)e2 

= r X mv H r x (r x s). (29) 

2mc r"^ 

The rate of change of L is given by 



■ r X (r X s) 



(.9-l)e 



2 



= I + [vx{r X s)+rx(r x s)+rx{v x s) — (3/r ){r-v)rx{r x s)], (30) 
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where I is the rate of change (|27p of the kinetic angular momentum r x mv. Using 
(I2Z1) in dSOl) gives 

• (g~ l)e^ r / X / M 

L = -— — -- [r X [V X s) + V X [s X r)\ 

Zmc^r-^ 

= -'-l-irir--)-s. (31) 

Combining this result with the Thomas rate of change of spin (|23p yields 

L + s = 0. (32) 
The total angular momentum L + s is thus conserved. 



4. Concluding remarks 

It is interesting to note that Thomas deduced in his landmark paper |2J that the orbit- 
averaged value of the angular momentum i + s of a spinning electron is conserved in 
its motion in a Coulomb field. Since Thomas's analysis, little attention seems to have 
been paid to the problem of angular momentum conservation in semiclassical models 
of hydrogen-likc atoms. Recently, however. Lush has revisited this problem and, using 
a different approach to ours, concluded that not even the orbit-averaged value of Z + s 
is conserved in this problem j25l [26]. Lush attributes the reason for the different 
finding in [5] to Thomas's omission of the contribution of hidden momentum to the 
force on a magnetic dipole and suggests the Thomas precession as the source of the 
nonconservation. 

In the problem at hand, omitting the Thomas precession term in a nonrelativistic 
Lagrangian would result in the replacement of the factor g — lin the spin-orbit coupling 
strength and the rate of change of spin by just the factor g. This, for the electron's 
g-factor of 2, would result in a spin-orbit strength by a factor of 2 greater than the 
empirically correct one and a rate of change of spin twice as great as that of the 
nonrelativistic limit of the Thomas equation. Moreover, irrespective of these fiaws, 
the angular momentum Z + s of the orbiting electron would still not be conserved in 
general. 

We note that laws of conservation of angular momentum hold true in general only 
when the angular momentum of the electromagnetic field in the system concerned is 
included [27l[28]. However, it can be shown easily that in the problem considered here 
the angular momentum of the central Coulomb field and the magnetic field of the 
orbiting particle's intrinsic magnetic moment, defined with respect to the Coulomb 
field's centre, vanishes in the nonrelativistic limiQ The angular momentum of the 
central Coulomb field and the magnetic field due to the particle's orbital magnetic 
moment can be shown to vanish, too [30] . 

The finding of our analysis is that the orbital quantity that is relevant to angular 
momentum conservation in a classical system with spin-orbit interaction is an orbital 
angular momentum defined in terms of the canonical momentum, not the standard 
kinetic momentum. Only the sum of such orbital angular momentum and the spin is in 
general conserved in the motion of a classical charged particle with spin in a Coulomb 

+ The angular momentum of the fields of a magnetic dipole /x and a charge q with respect to the 
position of the dipole is given by Mf = qr X {fj, X r)/cr^ , where r is the charge's displacement from 
the dipole 29 . But the field angular momentum M'^ about the position of the charge vanishes since 
M'^ = Mf — r X P{, where Pf = q{fi X r)/cr^ is the field linear momentum |29| . See also |27| . 
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plus the associated spin-orbit potential. This result is a forceful reminder of the fact 
that the quantum-mechanical operator of momentum corresponds to the canonical 
momentum of classical mechanics. This fact alone ensures that there is harmony 
between the quantum-mechanical and classical treatments in such an important aspect 
of the problem considered here as the conservation of angular momentum. 
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Appendix A. 

Here, we show that force (|16p with E = er/r^, 

F = -(fj, •V) (-V X e] +-E X 

\c / c at 

e , r e dn r , . . 

= --x(M-V)^---x-, (A.1) 

equals the spin-dependent part of force (fT3|) with Vso{r) = ge^ /2m'^c^r^ and s — 
-(2TOc/ge)/x, 

Fs = —■m '^^^°^^^ [r • (v X s)]r — mVso(r)(v x s) + — [mVso(r)(s x r)] 
ar at 

3e r / M - e , , e , „s r e dii r / . „n 

= 4 k • X fi)]r + —{v X ti) - - fi X {v V)— - - -77 X — . (A.2) 

cr cr"^ c r^ c at r 

To this end, we use the differentiation formula 

("•^)J = ^-3(a.r)J (A.3) 

to perform the following evaluation: 

r r 2 3 
V X {fi • V)— - fi X {v • V)— ^ —{v X fi) g-[(At • r)v - [v • r)fi] x r 

2 3 

— {v X fi) -[r X {v X fi)] X r 

1 3 
— -{v X fi) + — [r • {v X fi)]r. (A.4) 



Using this result in (jA.ip yields (|A.2|) . 



Appendix B. 

Here, we show that a symmetry of Lagrangian ^ implies the conservation of the total 
angular momentum defined as the sum of the canonical orbital angular momentum 
and spin. 

The change 6C of this Lagrangian under the rotation of the orbital and spin 
degrees of freedom qi = r , q2 ~ a, = b and their time derivatives about the origin 
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r = by an infinitesimal oriented angle 6<p, so that 6qi = Scf) x Qi and 6qi — 6(f) x q^, 
i = 1,2, 3, is given by 




3 

= ^[Pi- i^4> X Qj) + Pi • (<50 X Qi)] 

i=l 

3 

= 54)-^{q^xp, + q,xpi). (B.l) 

Here, in the second line, we introduced the canonical momenta pi = dC/dqi and used 
the Lagrange equations pi = dC/dqi, i = 1,2,3. Since Lagrangian ([3]) is manifestly 
invariant under a rotation of the vectors r, a, b and their time derivatives, we must 
have that SC = 0, and thus 

3 d ^ d 

= XI ^ Pi + Qi X Pi) ^ -^Yl 9' '^^'^ di^^^ ^^''^^ 

i—l i—1 

where we used the facts that qi x pi = L, which is the canonical orbital angular 
momentum (P^)) . and q2 x p2 + qs x p:^ = a x b = s, which follows from p7|) and ^TE\\ . 
This proves that the angular momentum L + s is conserved. 
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